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This paper investigates the parameter space of theories with gauge mediated supersymmetry 
breaking leading to gravitino (cold) dark matter with mass m 3 / 2 ~ IkeV —> lOMeV. We pay 
particular attention to the cosmological role of messenger fields. Cosmology requires that these 
messengers decay to the visible sector if the lightest messenger mass Mx 30 TeV. We then 
examine the various possible messenger number violating interactions allowed by the symmetries of 
the theory and by phenomenology. Late messenger decay generally results in entropy production 
hence in the dilution of pre-existing gravitinos. We find that in 5(7(5) grand unification only specific 
messenger- matter couplings allow to produce the required amount of gravitino dark matter particles. 
Gravitino dark matter with the correct abundance is however expected in larger gauge groups such 
as 5*0(10) for generic non-renormalizable messenger- matter interactions and for arbritrarily high 
post-inflationary reheating temperatures. 

PACS numbers: 98.80.Cq, 04.25.Nx 



I. INTRODUCTION 

Supersymmetric extensions of the standard elec- 
troweak model of particle physics come with a variety 
of appealing byproducts, such as the stabilization of the 
Higgs mass against radiative corrections, radiatively in- 
duced electro-weak symmetry breaking at the electro- 
weak scale, and the possibility of SU(3) C X SU(2)l x 
U(1)y gauge coupling unification at a sufficiently high 
energy scale. Most of these features occur rather nat- 
urally after supersymmetry (SUSY) has been broken 
softly, and it is believed that ultimately such a breaking 
has to occur spontaneously within some theory describ- 
ing all four fundamental interactions. 

Despite the lack of a particularly compelling model, 
there is a number of proposals with various theoreti- 
cal and phenomenological merits, where the spontaneous 
SUSY breaking usually takes place dynamically in a (hid- 
den) sector of the theory which does not contain the stan- 
dard model particles. Models can be classified according 
to the origin of SUSY breaking and of the soft terms, 
i.e. how the breaking of supersymmetry is transmitted 
to the low energy (visible) sector. In the so-called gauge 
mediated supersymmetry breaking (GMSB) models, this 
transmission is induced by renormalizable gauge interac- 
tions 0] (see Q for a review) . Particularly attractive 
features of these scenarios are the natural suppression of 
neutral current flavor changing interactions as well as a 
highly predictive mass spectrum that will be put to test 
in forthcoming collider experiments. 



In theories with gauge mediation, the lightest super- 
symmetric particle (LSP) is the gravitino and its mass 
can lie anywhere in the range 7/13/2 ~ 1 eV — > 1 GeV. 
Such a light gravitino (i.e. when TO3/2 1 keV) is tra- 
ditionally associated with a cosmological catastrophe. 
Indeed many studies have examined the production of 
light gravitinos in the early Universe in order to place 
stringent bounds on the post-inflationary reheating tem- 
perature (see 0, [H Q and references therein). Few 
studies have contemplated the possibility that this grav- 
itino LSP could make up the dark matter of the Uni- 
verse 

S U G3 IH E2, mainly because in the most 
naive model one needs to adjust the reheating temper- 
ature as a function of the gravitino mass in order to 
obtain the required dark matter relic density. However 
Refs. 0, have recognized the important cosmological 
role of the messenger particles that are part of the spec- 
trum of all GMSB theories. In particular, it has been 
shown that the late decay of the lightest messenger to 
visible sector particles can induce a substantial amount 
of entropy production which would result in the dilution 
of the predicted gravitino abundance. As a result, the 
light gravitino problem could be turned into a light grav- 
itino blessing, i.e. one would obtain suitable gravitino 
dark matter for arbitrarily high reheating temperatures. 

These studies Q, have focused on two specific cou- 
plings between the messenger and visible sectors. This 
motivates us to examine in more generality the possibility 
of producing the right amount of gravitino dark matter in 
GMSB scenarios. We do so in the present paper by con- 
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sidering all messenger-matter interactions allowed by the 
gauge symmetries of the theory and by phenomenology 
and by considering their impact on the gravitino abun- 
dance. The present study thus aims at being more ex- 
haustive than prior investigations; on the way we will also 
improve on some results previously obtained. In particu- 
lar we show that the coupling introduced in Ref. Q does 
not appear in minimal GMSB models and that multi- 
goldstino production channels modify substantially the 
results of Ref. @ . Finally we also take into account strin- 
gent constraints from big-bang nucleosynthesis and large 
scale structure formation. 

Our study is similar in spirit to those conducted for 
neutralino dark matter in minimal supergravity |13| . 
However since the present paper is of an exploratory na- 
ture, we approximate the mass spectrum of GMSB mod- 
els by three parameters: the messenger mass scale Mj, 
the supersymmetric particles mass scale Msusy ~ 1 TeV 
and the gravitino mass I n particular, we treat TO3/2 

and Mx as the fundamental parameters in our search for 
gravitino dark matter. We also discuss the influence of 
the nature and mass of the next-to-lightest supersym- 
metric particle (NLSP). Furthermore we calculate the 
velocity of the dark matter gravitinos in order to ex- 
amine whether this dark matter is hot, warm or cold. 
For simplicity we assume that R— parity holds. However, 
we will show that, when m 3 / 2 ^ lOMeV, our results re- 
main valid even if R— parity is violated. In this range the 
gravitino lifetime becomes much longer than the age of 
the Universe, so that it can be considered as stable on 
our cosmological time and constraints from diffuse back- 
ground distortions are eluded El E3. Finally we note 
that Ref. 0, 0] has discussed the possibility of grav- 
itino dark matter in minimal supergravity models (with 
gravity mediation of supersymmetry breaking). The cos- 
mology of these models is different from that of GMSB 
scenarios as the gravitino is heavier (m 3 / 2 > lOGeV) 
and there are no messenger particles. It is found in these 
studies that only a limited region of parameter space can 
satisfy the big-bang nucleosynthesis constraints and that 
the reheating temperature must be tuned in order to ob- 
tain the required dark matter relic density. 

The plan of this paper is as follows. In Section ITT1 we 
review the basics of gauge mediation models and discuss 
the nature of the lightest messenger particle which plays 
a crucial role in our analysis. In Section IlIII we discuss 
the numerous channels of light gravitino production in 
the early Universe as well as gravitino dilution due to 
late decay of a massive particle and the relevant cosmo- 
logical constraints. In Section llVI we survey the various 
renormalizable and non-renormalizable messenger num- 
ber violating interactions allowed by the gauge symmetry 
of the theory and discuss their consequences with respect 
to the light gravitino problem and gravitino dark matter. 
Finally in Section we discuss various perspectives, no- 
tably with respect to the case of SO(10) grand unifica- 
tion. We restrict ourselves to GMSB scenarios in the 
framework of N = 1 D = 4 supergravity and use natural 



units % = c = &b = 1; "ipi — 2.4 x 10 18 GeV denotes the 
reduced Planck mass. 



II. MESSENGER SECTOR 

A. Gauge mediation of SUSY breaking 

Gauge-mediated SUSY breaking (GMSB) 0,0], is usu- 
ally implemented by adding a term 

W = S$M$M + AW(S,Zi) (1) 

to the superpotential where <&m and $m are messen- 
ger left chiral superfields with SU(3) x SU(2) x U(l) 
quantum numbers, whereas the spurion left chiral super- 
field S and the secluded sector Z t fields are electroweak- 
and strong- interactions singlets. Upon the develop- 
ment of a non-vanishing vev (S) of the scalar compo- 
nent of the spurion superfield and a SUSY-breaking ex- 
pectation value of the spurion auxiliary field Fs, due to 
unspecified dynamics in the secluded sector AW(S, Z), 
fermionic messengers combine into Dirac fermions of 
mass M x ,\/2 — Mx = (S), whereas their bosonic part- 
ners mix in a mass matrix of eigenstates 4> and <\> with 
masses M = M X {1 - Fs/Mj,) 1 / 2 and A% = M x [l + 

Fg/M^-) 1 / 2 . In terms of the messengers bosonic compo- 
nents, = (-$* M + $m)/V2 and ~4> = ($ M + *m)/\/2. 
Note that 4> denotes a set of scalar fields transforming un- 
der some representation of the grand unified gauge group, 
and that the mass degeneracy of this multiplet is lifted 
by I?— terms and radiative corrections. 

Since the messengers share the standard model gauge 
interactions, the gaugino and scalar spartners acquire 
mass at the one- and two-loop levels respectively: 

hence the quantity A = F$/Mx is the supersymmetry 
breaking scale in the visible sector. Provided F$/Mx ~ 
100 TeV, this generates the required order of magnitude 
for the soft parameters 1 . Note that Fs < M x is manda- 
tory otherwise one of the messengers bosons acquires a 
negative mass squared. This also implies Mx ^ 100 TeV, 
and for M x > 100 TeV, F s < M' x , hence M x sets 
the mass scale for the messenger sector. In particular, 
w Mx, and in the following no distinction will be 
made between these two mass scales, except where oth- 
erwise noted. 

The gravitino mass is related to the fundamental SUSY 
breaking scale, to 3 / 2 = Fj \/3mpi, with F = Fs +Yli -^Zi 



In the present exploratory study, we do not address the detailed 
features of the GMSB mass spectrum and the related electroweak 
symmetry breaking and fine-tuning issues. 
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the sum of F— terms in the secluded sector. We define the 
parameter k = F s /F < 1, so that m 3 / 2 = Fg / (ky/3mp\) ■ 
In direct gauge mediated scenarios, one expects k < 1, 
whereas in scenarios in which the transmission of super- 
symmetry breaking to the messenger sector is loop sup- 
pressed one may find k <C 1. Note that one can also 
relate the parameters k, Mx and to 3 / 2 via the following 

formula: m 3 / 2 = AMx/ (kV^mpu- Since A is tied to 
the electroweak scale, the latter equation allows to elim- 
inate one parameter, which we choose to be k, in terms 
of TO3/2 and Mx, which we will treat as the fundamental 
parameters. 



B. Lightest messenger 

Taken at face value, GMSB scenarios generically lead 
to a cosmological catastrophe, as they predict that the 
lightest messenger should overdose the Universe 2 . In ef- 
fect, messenger gauge interactions as well as those de- 
rived from Eq. Q conserve messenger number so that the 
lightest messenger (a boson) is stable in this minimal ver- 
sion of the theory. As messengers can be easily produced 
in the primordial plasma thanks to their gauge interac- 
tions, their present day abundance is given by the result 
of a thermal freeze-out of messenger annihilation (akin 
to the well-known neutralino LSP freeze-out in gravity 
mediated SUSY breaking). 

Through explicit computation, one can show that the 
lightest messenger generically overdoses the Universe un- 
less its mass Mx ^ 10 4 GeV [l^. By lightest messen- 
ger, it is understood the lightest component of <fi af- 
ter taking into account D— terms and radiative correc- 
tions. Henceforth we denote this component by X. It 
has been shown that if the messengers sit in 5 + 5 rep- 
resentations 3 of SU(5), the lightest messenger carries 
the gauge charges of a sneutrino z>£, and its relic abun- 
dance would be of the right order of magnitude provided 
its mass - 10 - 30TeV Note that the mass 

scale in the messenger sector is a priori unconstrained, 
since phenomenology constrains the ratio Fg /Mx as dis- 
cussed previously. If the messengers sit in 10 + 10 repre- 
sentations of SU(5), the lightest messenger carries the 
gauge charges of a selectron Sr. Charged dark mat- 



2 Obviously, similar problems can in principle arise also for the 
Zi fields present in if the lightest Zi mass is larger than 
~ 100 TeV 1 181 . However, the issue becomes much more model- 
dependent here, and we will thus assume for simplicity that se- 
cluded sector fields can decay rapidly to the spurion field S. The 
latter field is free from such cosmological problems: even though 
it can be either heavier or lighter than the messengers [see the 
discussion following Eq. (3)] , in the first case it decays at tree- level 
to messengers, while in the second it decays to gauge bosons and 
gauginos fairly quickly through one-loop effects which are not 
suppressed by supersymmetry. 

3 Messengers sitting in complete GUT representations preserve au- 
tomatically gauge coupling unification. 



ter is forbidden by cosmology |20j and moreover this 
messenger would overdose the Universe for typical val- 
ues of Mx- Finally, if the messengers sit in 16 + 16 
representations of 5*0(10), the lightest messenger is a 
£ fi -like 517(3) x SU(2) x u ( l ) singlet. The next-to- 
lightest messenger can only decay to the lightest messen- 
ger by GUT scale suppressed interactions, and its life- 
time of order 10 10 yrs (M x /100TeV) -5 is so long that its 
decay produces unacceptable distortions of the diffuse 
backgrounds [l8j|. 

Cosmology thus forbids the lightest messenger to be 
stable unless messengers can be diluted to a very low 
abundance or the lightest messenger happens to have 
mass Mx ~ 10 — 30 TeV. If the post-inflationary reheat- 
ing temperature is larger than Mx and no late-time en- 
tropy production occurs, then messenger number must 
be violated, i.e. the Lagrangian of the theory must con- 
tain additional messenger-matter interactions. However 
such terms can spoil the phenomenological successes of 
the minimal model, in particular the absence of flavor 
changing neutral currents or an adequate pattern of elec- 
troweak symmetry breaking. One is thus tempted to be- 
lieve that such further messenger interactions with visible 
sector particles are rather weak, possibly resulting from 
non-renormalizable operators. This will be discussed in 
more detail below. 

Delayed messenger decay can have dramatic conse- 
quences for the gravitino problem and/or the possibil- 
ity of gravitino dark matter. If a non-relativistic species 
comes to dominate the energy density of the early Uni- 
verse and subsequently decays into visible sector parti- 
cles, a secondary epoch of reheating results and is con- 
commitant with the dilution of any pre-existing relics, 
such as gravitinos. The abundance of these relics may 
then, even for "arbritrarily" high primary reheat tem- 
peratures of the Universe after an inflationary epoch, be 
in accord with current observational constraints Q, H, • 

A crucial element in this analysis is the messenger 
abundance before decay. This is given, as mentioned 
earlier, by the thermal freeze-out of messengers, hence 
by their annihilation cross-section. The mass splitting 
in the messenger multiplet is generally small, of order 
Fs/Mx <C Mx, hence one should in principle consider 
the various co-annihilation channels. This task is how- 
ever left to a future more refined study; we note that over 
most of parameter space the inclusion of co-annihilation 
channels should modify the relic abundance of the light- 
est messenger by at most a factor of order unity since 
the the various particles that would co-annihilate have 
comparable annihilation cross-sections [l8[ . 



1. Annihilation cross-section 

Annihilation through gauge interactions. In the case 
of SU(5) unification, Dimopoulos et al. have calcu- 
lated the annihilation cross-section of the lightest messen- 
ger through gauge interactions for 5 + 5 representations 
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and parametrized it as: 



{(JXX'V) 



Ml 



x 



(3) 



with x = M x /T and A ~ —B ~ 3-10 -3 . This calculation 
takes into account all annihilation channels (mediated by 
gauge interactions) to gauge bosons, Higgses, neutrali- 
nos, charginos, fermions and sfermions, see Ref. [l8j| . 

Annihilation into two goldstinos. The above calcula- 
tion neglects annihilation into two goldstinos XX* — > 
GG. This latter occurs through a variety of diagrams: 
in the t— and u— channels, the annihilation takes place 
through the exchange of the fermionic mass eigenstate 
partner of X. In the s— channel, the annihilation oc- 
curs through the exchange of a graviton, a spurion and 
other secluded sector scalar particles. Finally, annihila- 
tion also occurs through four-point contact interactions 
XX*GG. These various contributions are triggered by 
various operators in the supergravity Lagrangian, taking 
into account the goldstino component of the gravitino 
and the fermionic spurion %ps fields after supersymmetry 
breaking as follows: 



/2 <9uG 



3 TO 3 / 2 



F s - 



(4) 
(5) 



For instance, the four-point contact interactions XX*GG 
derive from the gravitino mass term in the supergravity 
Lagrangian after expanding the exponential of the Kahler 
function to second order in the scalar around their vev. 
The Yukawa coupling between the goldstino and the hid- 
den sector SUSY breaking scalars Zi, which enters the 
s— channel exchange diagrams, also derives from this ex- 
pansion. The coupling XX*Zi derives from the scalar 
potential trilinear couplings. Finally, the coupling be- 
tween X, its fermions mass eigenstate partner and G is 
obtained directly from the supergravity Lagrangian cou- 
pling between the gravitino and a pair of fermion-boson 
partners, taking the appropriate linear combination to 
express it in terms of the mass eigenstates after SUSY 
breaking. 

The annihilation cross-section into two goldstinos must 
be calculated with care since some leading high energy 
contributions are expected to cancel out |2l], |2^, |2^| . It 
will be important to distinguish between the cases where 
the spurion mass Ms is larger or smaller than that of the 
lightest messenger. Both configurations are dynamically 
possible: for instance, one finds in the simplest models |2j 
that M\ = (k — \/3A)(k/A 2 )M|, where A and re denote 
respectively the spurion self-coupling and its coupling to 
the messenger fields in the superpotential, and where we 
have neglected here the effect of the spurion coupling to 
the secluded sector fields following the study of [24| for 
the stability and local minima conditions. In this case 



the spurion is heavier than the lightest messenger when 
\/3 < re/A < 2.2 and lighter when re/A > 2.2. 

In the light scalars and non-relativistic limit, s ~ 
4M X (1 + 3/x) 3> m%. with mz ; mass of the secluded 
sector scalar Zi, and x = Mx/T as above, the cross- 
section reads: 



'xx* 



Fs 
F 



1 



32ttM x 



In effect, performing an expansion in Fs/s in the matrix 
element, one finds that both terms of order and 1 cancel 
among the various contributions, yielding a cross-section 
ex Fg. Its high energy limit s ^> M x , m%. is actually that 
of spurion-spurion annihilation into two goldstinos |23j |. 
In practice, Eq. © applies if the spurion is much lighter 
than the lightest messenger since its coupling to XX* 
dominates that of secluded sector scalars as a result from 
the tree-level coupling in the superpotential between S, 
$m and 'I'm. However, in this limit one must also ac- 
count for annihilation of X, X* into spurions; this will be 
discussed further below. 

If secluded sector scalars are heavier than the lightest 
messenger, the limit s <C m|. applies and in this case, 
one can neglect the s— channel exchange of these scalars. 
The cancellation between the various diagrams occurs 
only to order in (Fs/s), leaving a cross-section oc Fj: 



'xx* 



*GG V ) 



1 FlM 2 x 
8tt F 4 



2x 



(7) 



The s— channel exchange graphs of the secluded sector 
scalars are suppressed by s/m 2 z .. Note that the term Fs 
in these expressions should be understood as the mass 
squared difference between the fermion and boson com- 
ponents of the lightest messenger multiplet, rather than 
as the vev of the auxiliary component of S. These two 
quantities differ if the superpotential includes a coupling 
constant re, W D kS$>$>; our choice here is re = 1. 

The annihilation cross-section in the heavy spurion 
limit [Eq. Q] violates the unitarity bound (cr ann u) < 
8tt/M x [251 (in the non-relativistic regime) for Mx ^ 
1.63 • 10 7 GeV (m 3/2 /l keV) 2 / 3 . Beyond this limit the ef- 
fective Lagrangian is no longer valid, and one expects 
sizeable contributions from multi-goldstinos production. 
Hence the results obtained hereafter in the region where 
the unitarity bound is violated are highly uncertain and 
model-dependent. In what follows, we assume that the 
cross-section saturates at the unitarity limit in this re- 
gion M x > 1.6 • 10 n GeV (m 3/2 /l GcV) 2 / 3 if the spurion 
is heavier than the lightest messenger. We also consider 
the other possible limit in which the cross-section follows 
Eq. ©, so that the comparison of these two cases will 
allow us to assess the impact ofthe above effects on the 
relic gravitino abundance. 

Since phenomenology requires that Mx ~ Fs /10 5 GeV 
(see Section II. A), it is easy to see that annihilation 
into goldstinos dominates the cross-section in the heavy 
scalars limit for M x > 3 • 10 6 GeV (m 3/2 /l keV) 2 / 3 . 
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Hence the inclusion of this channel in the present cal- 
culation modifies rather drastically the relic abundance 
of the lightest messenger in this part of parameter space. 

Annihilation into two spurions. If S is lighter than X, 
there is no problem associated with unitarity, and one can 
safely use Eq. © all throughout parameter space. One 
must nonetheless account for XX* — > SS annihilation, 
whose cross-section reads 



(o~xx'->ssv) 



1 



QAttM x 



1 



1 . 3 
x + ^2 



23, 
Ax' 



(8) 



to first order in x~ x and in r{= Alg/M x ), and where we 
neglected for simplicity the contribution of a AS 3 term 
in the superpotential, assuming that A -C k ~ 1. 

This cross-section is comparable to the annihilation 
cross-section through gauge interactions given in Eq. (J3J 
and results in the decrease of the relic abundance of the 
lightest messenger by a factor of order 2. For direct 
GMSB models in which F$ ~ F, the annihilation chan- 
nel into goldstinos becomes dominant and must be taken 
into account. 



2. Relic abundance 

Freeze-out of the lightest messenger annihilations oc- 
curs at a value x / : 



x f ~ log 



0/1 



B/A \ 1 

i°g(Q/)y v^gioj) 



(9) 



where Q f ~ 6.1 • 10 10 (M X /10 6 GeV)- 1 ^, and the values 
of A and B accounts for the various possible channels 
depicted above. In terms of this freeze-out value Xf, the 
relic abundance is then given by: 



Y x ^ 2.1 • 10 



-14 ( My 



V10 6 GeV / A + B/2x f ' 



In the case of messenger sitting in 10 + 10 represen- 
tations of SU(5), the relic abundance is expected to be 
similar to the above to within a factor of a few, since the 
lightest messenger carries hypercharge. For simplicity, 
we thus assume that its relic abundance is also given by 
Eq. JTJJ above. 

Finally, in the case of SO(10) grand unification, the 
lightest messenger is a singlet under the standard model 
gauge interactions. As argued in Ref. [26j, it can anni- 
hilate through one-loop diagrams (which dominate the 
exchange of tree level GUT mass bosons considered in 
Ref. [l8j) and into two goldstinos at tree level as above. 
This case has been discussed in some detail in the low 
Mx region in Ref. j2|| . In Section[V]we sketch briefly the 
parameter space of 50(10) GMSB scenarios using order 
of magnitude estimates of these diagrams. In the main 
discussion of this paper, we thus focus on SU(5) grand 



unification with messengers either in 5 + 5 or 10 + 10 
representations . 

The cosmological scenario we have in mind is the fol- 
lowing. As the Universe reheats to high temperature af- 
ter inflation, radiation along with gravitinos and messen- 
gers are produced. As the temperature decreases, the 
lightest messenger annihilations cease and its abundance 
freezes-out. This non-relativistic lightest messenger may 
come to dominate the energy density if its decay to the 
visible sector is sufficiently delayed and its relic abun- 
dance sufficiently large. In particular, messengers come 
to dominate the energy density when the background 
temperature 



T dom ~ -M X Y X 



(11) 



(io) n 



(provided Tdom ^ Ticc, with Td cc the temperature at 
which messengers decay, see below). During decay of the 
lightest messenger, gravitinos and possibly sparticles may 
be produced, but the pre-existing gravitino and NLSP 
abundances are diluted by entropy production. Finally, 
at late cosmological times, the NLSP decays to the grav- 
itino. In some cases the NLSP may decay before the 
lightest messenger. The final gravitino abundance is then 
the sum of gravitinos produced by sparticle and messen- 
ger interactions at early times and diluted later by the 
appropriate factor due to messenger decay, plus graviti- 
nos produced during the secondary reheating induced by 
messenger decay as well as gravitinos produced in the 
decay of the NLSP. Then, for a given messenger decay 
width, which characterizes the dilution factor, one may 
find in the to 3 / 2 — Mx plane the region where satisfactory 
abundances for gravitinos are found. 

If Tciom < Tdcc, the lightest messenger never comes 
to dominate the energy density and no entropy produc- 
tion ensues. This is what has been generally assumed in 
previous studies that derived upper limits on the post- 
inflationary reheating temperature from the upper limit 



3/ 2 < 1, e.g. 
tion IV.A.l. 



This case is discussed in Sec- 



III. GRAVITINO PRODUCTION 

A. Production channels 

The fractional contribution to the present critical den- 
sity of non-relativistic gravitinos G with number-to- 
entropy ratio Y3/2 = n/s is given by 



f> 3/2 fr 2 = 2.81-10 8 (y 



m 3 / 2 



GeV 



Y ? 



3/2 



(12) 



where h is the Hubble constant in units of 100 
kms~ 1 Mpc~ 1 . The number-to-entropy ratio Y 3 / 2 is 
found by following the Boltzmann equation describing 
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gravitino production: 



dY 
dT 



sHT ^ 



) +22^ a i+j^G+... V V n i n j) 



(13) 

which includes production by sparticle decays , 
and scatterings (a i+ ^Q + VijUiUj), neglecting three- 
body and higher order interactions. In principle, one 
should take into account gravitino annihilation as well. 
However it is sufficient to approximate possible grav- 
itino losses by imposing that the gravitino number-to- 
cntropy ratio never exceeds its thermal equilibrium value 
Y cq ~ 3.7 x 10~ 3 (g*/230)" 1 , with g± the number of rela- 
tivistic degrees of freedom in the thermal bath 4 . 

The production of gravitinos in the early Universe 
is dominated by the production of the helicity ±1/2 
component if m 3 / 2 *C M, where M denotes the mass 
scale of particles leading to gravitino production. This 
helicity ±1/2 component is related to the goldstino 
through Eq. J2J, associated with the breaking of local 
SUSY, where ^ and G denote respectively the helic- 
ity ±1/2 gravitino components and the goldstino |27j . 
In scenarios of gauge mediated SUSY breaking in which 
™3/2 ^ 1 GeV<C Mstjsy ~ ITeV, the correspondence 
^ ~ iy/2/Z 9 M G/m 3 / 2 , is generically satisfied. The 
production of goldstinos may then be decomposed into 
the contributions from particle scatterings, decay of par- 
ticles before the freeze-out of the NLSP, and the possible 
contribution of NLSP decay. 

The decay and scattering contributions of visible sector 
fields have been calculated in various studies using the 
above gravitino-goldstino equivalence (see e.g., Ref. Q 
for a detailed discussion and references). Later on it has 
been shown that the contribution of messengers is quite 
significant [Sj- Indeed messengers couple directly to the 
goldstino G through the superpotential term S$m$M 
and the fraction of goldstino G comprised in the fermionic 
component of S, Eq. (JSJl. We use the decay widths and 
cross-sections for messengers interactions leading to grav- 
itino production given in Ref. 0. 

These authors have also argued that the goldstino de- 
couples from the visible sector fields at energies Mx ^ 
E < Vf claiming that the effective gravitino-particle- 
sparticle vertex is induced by loop diagrams involving 
messengers |28|. However this analysis has been per- 
formed in the limit of global supersymmetry and it 
thus ignores the tree level fermion-sfermion-goldstino 
- (1 / \/2mp\)gij* d f _ t <&* 3 XRl u J tJ '^ ' v +h..c. and gauge boson- 



gaugino-goldstino (i / Am Pi )^ ^ ^ [F$ +F f 



'A a h 



pa 



interaction terms which appear in local supersymmetry. 
Therefore the goldstino does not decouple from visible 



4 Note that this equilibrium abundance only includes spin 1/2 
goldstinos, with the population of 3/2 components of the grav- 
itino assumed to be negligible. 



sector fields at energies Mx ^5 E < \/F, at variance with 
Ref. @- Nevertheless, it is true, following Refs. [ft l2?H 
that production channels should include the loop-induced 
messenger contribution to the effective particle-sparticle- 
goldstino vertex at energies E < Mx- This is done in 
the present study; we do not include thermal corrections 
to the cross-sections since they are found to be negligi- 
ble [13. 

The various contributions to Y 3 / 2 have different tem- 
perature dependences. The largest fraction of the decay 
contributions results from cosmic epochs when the tem- 
perature falls below the mass of the corresponding parti- 
cles. In contrast the dimension-5 operators associated to 
visible sector particle - sparticle scatterings are most ef- 
fective at high temperatures; the contribution from mes- 
sengers peaks at temperatures T ~ Mx- A simple fit to 
the results of Ref. @, with the modifications according 
to the remarks above, gives an estimate of the amount of 
gravitinos produced by scatterings and decays of sparti- 
cles and messengers for a reheating temperature Trh\ 
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(14) 



The first two lines correspond to sparticle decays and 
particle-sparticle scatterings, respectively, while the last 
two correspond to interactions involving messenger fields. 
Hence the latter should only be included if Trh ^ Mx- 
In these equations, M 3 denotes the gluino mass scale, 
which controls the amount of gravitinos produced. In- 
deed, the coupling between gauginos, gauge bosons and 
golstinos dominates the particle-sparticle-goldstino con- 
tributions at high temperatures and scales as the gauge 
coupling constant; its contribution is represented by the 
second term on the r.h.s. of Eq. 114[) . Decay contribu- 
tions to £l 3 / 2 h 2 scale as M 3 , with M the mass of the 
decaying sparticle, see the first term on the r.h.s. of 
Eq. I|14jl. The gluino appears in this term due to the 
larger degree of freedom for colored particles than for 
color singlets as well as due to colored sparticles being 
heavier than color singlets in the ratio a 3 jai or a 3 /a\ 
[see Eq. Strictly speaking this relation applies at 

the messenger scale and must be corrected by renormal- 
ization group running at lower energy scales. Neverthe- 
less the mass spectrum of GMSB models is dominated 
by squarks and gluinos whose masses are comparable at 
scale Mx- As regards the contribution of messengers, 
given in Eq. (|rl|) for Trh ^ Mx, the gluino mass ap- 
pears only in the combination (47r/a 3 )M 3 ~ Fs/Mx, 
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i.e. as the normalization of F$ as a function of Mx and 
the electroweak scale. 

We provide Eq. 114JI in order to assist the reader in 
interpreting the figures that follow. The results shown 
below are obtained from the integration of the Boltz- 
mann equation. Also recall that Eq. (|14J) does not take 
into account all contributions to the gravitino abundance. 
One must notably add the contribution of gravitinos pro- 
duced by annihilations of the lightest messenger, as well 
as lightest messenger decay and NLSP decay. 

The NLSP decays into final states including one grav- 
itino with width: 



JU NLSP 



NLSP^G 



48-7T m^i^m 



(15) 



temperature 



at 



NLSP^G 

-l 



The background 
NLSP decay is then T : 
5MeV(Af N LSp/100GeV) 5 / 2 (m 3/2 /lMeV)- i . This 
decay occurs late and consequently big-bang nucleosyn- 
thesis constraints on hadronic or electromagnetic energy 
injection at time ~ 10~ 2 — 10 8 sec lead to the exclusion 
of a significant part of parameter space, as will be seen in 
the following. Since one NLSP produces one gravitino, 
the gravitino yield of NLSP decays in terms of entropy 
density is simply Ynlsp- 

In principle, the NLSPs result from a freeze-out of 
thermal equilibrium, hence Ynlsp can be calculated in 
the same way as the relic abundance of the lightest mes- 
senger. However, one must not forget the possible pro- 
duction of NLSPs during the decay of the lightest mes- 
senger. If the decay temperature of the latter is larger 
than the NLSP freeze-out temperature, this contribution 
is washed out by NLSP annihilations. However if de- 
cay occurs later, NLSPs are regenerated to a level which 
depends on the time of messenger decay and on the anni- 
hilation cross-section of the lightest messenger. In more 
detail, if the lightest messenger decays to NLSPs before 
NLSPs decay in turn to gravitinos, the NLSPs have time 
to annihilate. In Ref. a procedure was outlined to 
calculate the number density of NLSPs remaining after 
these further annihilations. If the lightest messenger de- 
cays to NLSPs after pre-existing NLSPs have decayed to 
gravitinos, the calculation of the final number of graviti- 
nos produced is more involved and necessitates the inte- 
gration of coupled Boltzmann equations. For simplicity, 
we assume that all NLSPs produced in messenger decay 
then decay instantaneously to gravitinos without anni- 
hilating. This maximizes the number of gravitinos pro- 
duced hence reinforces the constraints from relic density 
arguments, and, in this sense, this assumption is conser- 
vative. Keeping track of NLSP annihilation before decay 
to gravitinos would not affect our results significantly as 
it is marginal in most of parameter space [8|. 

Finally there exist other potential channels of grav- 
itino production. One is that of helicity ±3/2 produc- 
tion, which for large gravitino mass and high reheat- 
ing temperature may become important. The helicity 
±3/2 modes interact with gravitational strength only 



and are produced by interactions in the thermal bath in 
abundance O^/V ~ (m 3/2 /lGeV)(r RH /10 14 GeV) 0. 
Therefore this contribution does not dominate in most 
of the m 3 / 2 — Mx — Trh parameter space. We do 
not take into account possible non-thermal production 
channels of helicity ±3/2 gravitinos during inflation [3l| . 
The amount of gravitinos produced in this way depends 
strongly on the underlying model of inflation, eventhough 
it may exceed the amount of helicity ±3/2 modes pro- 
duced by scatterings in the thermal bath in particular 
models. 



B. Gravitino Dilution 

The delayed decay of non-relativistic messengers may 
have dramatic consequences on the abundance of any pre- 
existing species, such as gravitinos. In case delayed decay 
results in the temporary matter-domination by messen- 
ger rest mass, i.e. px 3> pi where p denote energy densi- 
ties and subscripts X and "r" refer to messenger and radi- 
ation, respectively, entropy production is significant and 
results in the severe dilution of any pre-existing number- 
to-entropy ratio Y. In this case the post- messenger-decay 
cosmic radiation temperature T^ ec is substantially larger 
than the pre-decay temperature Tf~ , akin of a second 
reheat. Approximating decay to be instantaneous when 
the Hubble scale equals the decay width of the lightest 
messenger, H ss Tx, one finds 



J dcc 



xmpi, 



(16) 



where <?> denotes the number of relativistic d.o.f. at 
temperature . If the particle decays into the visible 
and into an invisible sector, the decay width in Eq. i|16|) 
above should be multiplied by -Bvisibie, with B V isibic the 
branching ratio into visible sector particles. This may be 
of relevance notably when the lightest messenger decays 
into visible sector particles and into gravitinos which do 
not share their energy density with the visible sector af- 
terwards. 

By equating the pre- and post-decay energy densities, 
the pre-decay radiation temperature is obtained in 
terms of T^ ec and Xd om [see Ea. lfTT)) ] at which X comes 
to dominate the energy density, as: 



J dcc 



T. 



dec 



9< 



1/3 



± dcc 
^dom 



1/3- 



(17) 



J dec 



, while if T> c « T d 

Olll * 



Obviously, if X does not dominate the energy den 
sity before decaying, T^ cc ~ T ,, ml. ;r r 
one finds T< c « T> 
substantial. In effect the ratio of pre-decay and post 
decay entropy densities, gives the entropy release Ax = 
s>/s < =g > Tl/g < Tl: 



dec ^dec an d entropy production is very 
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where it is understood that if Ax ^ 1 is found, it ought 
to be substituted by Ax = 1. 

If there exists a whole tower of N unstable, but long- 
lived particles, with abundances Yi, masses Mi and decay 
widths Tj for particle i, the final dilution factor is deter- 
mined solely by the properties of the slowest decaying 
messenger. In particular, all the equations above may be 
employed as if any prior decays had not occurred. 



C. Cosmological constraints 

1, Hot, warm or cold dark matter 



The values of and T^ om are given in terms of Yx, 
M x and T x through Eqs. CSJ and JTTJ. 

Such entropy release dilutes pre-existing densities ac- 
cording to: Y> = Y</Ax- Nevertheless, it should be 
borne in mind that, in case of high second reheat tem- 
peratures T^ oc , the regeneration of diluted species may 
occur. This effect is taken into account in our calcula- 
tions by treating messenger decay as a second reheat. 
Note that substantial entropy release after BBN is unac- 
ceptable, and > 1 MeV is required. Eq. i|16|l may 
thus be employed to infer a fairly strict lower limit of 
Tx ^ 4.3 x 10~ 25 GeV on abundant and slowly decaying 
particle species in the early Universe. 

Particularly interesting to cosmology is the case of 
significant entropy dilution. In this limit one may use 
Eqs. (|T%)) and ((T7|) to derive the entropy dilution factor 
for an abritrary species with mass Mx , decay width Tx 
and abundance Yx, in the limit Ax 1: 



One finds m > 1 keV from the req uirement that the Uni- 
verse has reionized by z ~ 6 [32( or from the measure- 
ment of the power spectrum in the Lyman a forest |33| . 
If reionization has occurred as early as z ~ 17, as sug- 
gested by the recent WMAP data [34J , then a mass larger 
than ~ lOkeV seems required [35) . 

These constraints are important to our analysis and 
we keep track of the average velocity extrapolated to 
2 = 0. Obviously the limits between hot, warm and 
cold matter are fuzzy, and we choose to qualify as warm 
dark matter particles with velocity 0.0018 kms -1 < 
vq < 0.054kms~ 1 , corresponding to particle masses 
0.3 keV < m < lOkeV (for freeze-out from thermal equi- 
librium as above). For velocities above the upper limit or 
below the lower limit, we mean hot or cold dark matter 
respectively. It is important to note that entropy produc- 
tion after decoupling of the gravitinos cools down these 

1/3 

dark matter particles according to: (vo) —> (vq)/A^- . 

For gravitinos produced by out-of-equilibrium pro- 
cesses, notably by the late decay of a massive particle, 
the above relation between mass and velocity is modifed. 
Assuming the outgoing gravitino carries a momentum of 
half the mass M of the decaying particle, the present 
velocity reads: 

(vq) « (Af/2m 3/2 )(3.91/ 5 *, doc ) 1/3 (T /T doc ), (21) 

1/4 

9*, doc the number of d.o.f. at decay, To the present 
cosmic background temperature and Td GC the tempera- 
ture at decay. Note that Td oc generally depends on M 
so that the dependence between the nature of gravitino 
dark matter (cold/warm/hot) and the mass of the decay- 
ing particle is not necessarily trivial. For instance, one 
can show that a decaying NLSP produces hot /warm dark 
matter if its mass < 500 GeV. 

If the decay occurs at temperatures sufficiently high 
that the gravitino can interact and thermalize, one should 
rather use Eq. 120|) . However at temperatures T <C 
100 GeV the gravitino has decoupled from the thermal 
plasma [2^|, mainly because the sparticles have decou- 
pled themselves. Therefore the decay of the NLSP or 
of the lightest messenger (if sufficiently late) generally 
produces highly relativistic gravitinos. 



Collisionless damping during the radiation era leads 
to the erasure of power in density fluctuations below a 
length scale (the free-streaming scale) that is mostly de- 
termined by the time at which dark matter particles be- 
come non-relativistic, or equivalently by their velocity 
extrapolated to zero redshift. A particle of mass m that 
thermally decouples from the plasma when relativistic 
has a present-day velocity: 

(vo) « O.OlSkms- 1 ( 5 ^ dcc /230)- 1/3 (m/lkeV)- 1 , (20) 

with g^dec the number of d.o.f. at decoupling. Cosmolog- 
ical data on the power spectrum of density fluctuations 
allow to place constraints on the mass of the particle. 



2. Big-bang nucleosynthesis constraints 

Due to the different channels of gravitino produc- 
tion, one generally finds that gravitinos are made of two 
generic sub-populations: one that has been produced by 
equilibrium processes and another made of hot graviti- 
nos produced by out-of-equilibrium decays. It may be 
that the latter are so highly relativistic that they form 
a hot dark matter component. However their impact 
on the formation of large scale structure may be negli- 
gible if their contribution to the gravitino energy den- 
sity is negligible. In this particular case, constraints 
from big-bang nucleosynthesis on extra degrees of free- 
dom may apply and constrain this population. We take 
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the BBN constraints on extra degrees of freedom to 
be Sg < 1.8 (corresponding to 1 extra neutrino fam- 
ily allowed) [2(j. Gravitinos that are relativistic at the 
time of BBN and carry energy density py 2 contribute 

to the level #<7bbn/<7bbn — P3/2I with p s r the stan- 
dard radiation energy density at the onset of BBN with 
5bbn = 10-75. Gravitinos that were once in thermal equi- 
librium or that were produced by equilibrium processes at 
temperature T (d.o.f. g) carry characteristic momentum 

P3/2 ~ 3T B bn(5bbn/5 s ) 1/3 at BBN, with T BBN ~ 1 MeV. 
Assuming g s ~ 230, these gravitinos are relativistic if 
m 3/2 ~ 1 MeV and their contribution to the energy den- 
sity is (5(?BBN /<?BBN — Y 3/2 witn Y 3/2 = n 3/2/ S ' nence ^ 

is negligible due to the upper bound on Y3/2 resulting 
from thermal equilibrium. 

However most relativistic gravitinos at the time of 
BBN result from out-of-equilibrium decays, e.g. from 
NLSP or from the lightest messenger decay. Given that, 
immediately after decay the outgoing gravitinos carry 
a fraction -B3/2 of the rest mass energy of the decay- 
ing particle of mass M as kinetic energy, their con- 
tribution to the energy density at the onset of BBN 
reads: <%sbn/<?bbn - S 3/2 My M (4/3)( 5 | jBN /< ? s ) 1 /3/T d , 
where Td is the decay temperature and Ym the number- 
to-entropy ratio of the parent at decay. Here it has 
been assumed that Td ^ 4MYjm/3. In the opposite 
limit, i.e. in the case of significant entropy production 
at decay, the above relation becomes <5c7bbn/.9bbn — 
-B3/2(Sbbn/5 s ) 1/3 (assuming B 3/2 < 1), and it will be 
this limit which results in the strongest BBN constraints. 

The time of decay of NLSPs to gravitinos is also 
strongly constrained by big-bang nucleosynthesis lim- 
its on hadronic and electromagnetic energy injection at 
times ~ 10~ 2 — > 10 8 sec. These constraints have been ex- 
amined in Refs. [II IH IH HI , while Ref - has trans- 
lated these bounds on the messenger scale Mx of GMSB 
models assuming k = 1 (which is equivalent to setting an 
upper bound on 7713/2)- in the present analysis, we use 
the latest constraints from hadronic and electromagnetic 
energy injection from Ref. |38l |. 

In GMSB scenarios, the NLSP is generically a neu- 
tralino (mainly bino) or a stau. The former decays pre- 
dominantly into a photon and a goldstino; the fraction of 
energy spent with branching ratio B cnl ~ 1 and by three 
body decays into a pair of quarks and goldstino with 
hadronic branching ratio Shad ~ 10~ 3 ; if its decay to Z 
bosons is not suppressed by phase space, i.e. (Mnlsp — 
Mz)/Mz ^ 1, the hadronic branching ratio -Bhad ~ 0.15. 
For simplicity, we use this latter value, which is conserva- 
tive with respect to the constraints inferred. Concerning 
the annihilation cross-section of the bino NLSP, we use 
the value (ctnlsp^) = lO^GeV" 2 (A/ N lsp/100 GeV)~ 2 , 
which corresponds to the bulk region of minimal super- 
gravity 0| . We will comment on the dependence of our 
results on the choices made when discussing the results 
shown in Fig.^below. 

A stau NLSP may produce in its decay electromag- 
netic and hadronic showers. About 100% of the energy 



is converted to electromagnetically interacting particles. 
In 70% of all decays, a stau NLSP produces hadrons, 
but these are mesons whose lifetimes are so short that 
they do not have time to interact before decaying if 
they were emitted at times > 10 2 sec. Hence, we use 
-Shad = 0.7 for stau decay timescales shorter than 10 2 sec 
and Shad = 10~ 3 for longer decay times. The stau anni- 
hilation cross-section is not as model dependent as that of 
the bino, (ctnlspw) - 10~ 7 GeV" 2 (M NLS p/100 GeV)~ 2 . 
Given its large annihilation cross-section, the stau has a 
small relic abundance, and consequently the BBN con- 
straints are comparatively weaker. 

Overall big-bang nucleosynthesis constraints apply to 
a combination of Inlsp (relic abundance) and decay 
timescale oc -^nlsp to 3/2- Note that for a decay timescale 
t ~ 10 3 sec, interesting modifications to BBN may re- 
sult [13. We also note that in a very limited part of 
parameter space of GMSB theories, the NLSP can be a 
sneutrino for which the BBN constraints would be 
largely reduced 0, . 

Finally, since the mass of the NLSP enters the BBN 
constraints while the gravitino yield is controled by the 
gluino mass, it is necessary to schematize the mass spec- 
trum of GMSB scenarios. We do so by assuming a 
mass ratio M 3 /Mnlsp ~ 6 d an d fiducial values 
M NLS p = 150 GeV and M 3 = ITeV. Where relevant we 
mention the possible influence of these values on our re- 
sults. 



IV. MESSENGER COUPLINGS TO MATTER 
AND GRAVITINO DARK MATTER 

In this section, we investigate the possible solutions 
for gravitino dark matter for various messenger number 
violating interactions added to the Lagrangian. As ar- 
gued in Section m such interactions are mandatory if no 
substantial entropy production occurs at temperatures 
< Mx (other than due to lightest messenger domina- 
tion and decay) in order to avoid the cosmological prob- 
lems that would result from the stability of the lightest 
messenger. For definiteness, we adopt the notations of 
Ref. [23 including four component spinors, in the gen- 
eral supergravity Lagrangian. 



A. Renormalizable couplings 

1. Superpotential couplings 

Renormalizable couplings, beyond those of the re- 
quired messenger gauge interactions, may exist, though 
they are constrained by considerations of flavor changing 
neutral currents as well as the potential development of 
charge- and color- breaking minima, among other issues. 
Dine et al. ^3 have analysed viable extensions of the 
minimal GMSB scenario in this direction, introducing 
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couplings of the form 

W D yiH D ^ l M e t + y\H D Q^ q M , (22) 

where & M and <& q M denote lepton- and quark- like mes- 
sengers. Here denotes an SU(2) messenger doublet, 
the 2/i's are Yukawa couplings with family index i, and 
Hd, e, and Qi are standard model down- type Higgs, 
right-handed lepton, and quark doublet, respectively. 

Additional SUSY-breaking mass splittings are gener- 
ated by these types of interactions via one-loop contribu- 
tions yielding, for example, negative relative mass con- 
tributions of order Smg/nii w — 10 3 yfFg/M x to slepton 
masses. Flavor changing neutral currents may place po- 
tentially restrictive limits on such couplings, as due to 
the experimentally verified weakness of such processes, 
mass splittings between 1st- and 2nd- generation slep- 
tons (and squarks) are constrained to be smaller than 
(m| — m| 2 )/?7i| < 10~ 3 . Assuming conservatively, 
yf 7^ and yf = yf = one may thus infer a limit 
y\ < (Mjsf/10 8 GeV) on this extra Yukawa coupling. 

Interactions induced by the superpotential Eq. lt2*2^) 
also induce the decay of messengers, in particular X — > 
H~e + assuming X carries the same gauge charges as a vl 
(see Section[nJ|. Hence the decay width Tx — y 2 Mx/8Tr. 
Though the limit on Yukawa couplings as inferred above 
may be quite severe, entropy production due to delayed 
messenger decay is absent when terms of Eq. (I22|l are 
included into W, unless the extra Yukawa coupling is 
extremely small, y < 10 _15 (M X /10 7 GeV) 3 / 2 , assuming 
SU(5) grand unification. 

It is nevertheless instructive to study the influence of 
such "fast" decay of the lightest messenger on the pos- 
siblity of having gravitino dark matter. As mentioned 
earlier, this case (without entropy production) has been 
implicitly assumed in previous studies that have drawn 
upper bounds on the post-inflationary reheating temper- 
ature from the upper bound on the gravitino density 
^3/2 < 1 0-E1 (and references therein). In order to pro- 
vide a point of comparison with this previous litterature, 
we plot in Fig. ^ the results of the calculation of f2 3 / 2 
in the plane Trh — using the techniques developed 

in the previous section. We assume "fast" decay with 
width Tx ~ 10~ 9 Mx corresponding to y ~ 10~ 4 , which 
ensures that phenomenological constraints are satisfied 
for all values of Mx- The results shown are insensitive 
to the exact value of y, provided it is not so tiny that 
substantial entropy production would occur. 

The shade (color) coding in this figure and all subse- 
quent figures is as follows: lightest (yellow) corresponds 
to f^3/2 < 0.01 (no gravitino problem but no dark mat- 
ter), and the increasingly darker (respectively green, red 
and blue) areas indicate respectively the regions of cold, 
warm and hot dark matter in which 0.01 < ^3/2 < 1- 
The area shaded by lines oriented NE-SW at the right of 
each figure corresponds to the region excluded by BBN 
constraints on NLSP to gravitino decay. White color 
indicates ^3/2 > 1, i.e., overclosure of the Universe by 



gravitinos. Finally the area marked with horizontal lines 
is unphysical as it corresponds to F$ > F (k > 1). 

We choose 0.01 and 1 as lower and upper bounds re- 
spectively to delimit where the gravitino can account for 
dark matter, eventhough cosmological data restrict this 
to a much smaller range. However the calculations pre- 
sented here contain intrinsic uncertainties of factors of 
a few that were mentioned in the previous sections, and 
therefore the green, red and blue areas should be under- 
stood as indicative of the region in which one can find 
solutions for gravitino dark matter. 

As indicated in the caption of Fig. ^ the left pan- 
els correspond to Mx — 10 5 GeV while the right pan- 
els correspond to Mx = 10 10 GeV (for which the con- 
dition k < 1 translates in m 3 / 2 250 keV). The upper 
and lower panels correspond respectively to a stau and a 
bino NLSP. As anticipated in the previous Section, the 
BBN constraints on hadronic and electromagnetic en- 
ergy injection do not apply to the stau at these small 
values of m.3/2, as a result of the low stau relic abun- 
dance. In fact, for a stau of mass Mnlsp — 150 GeV, 
a gravitino as heavy as m 3 ^ 2 ~ 10 GeV is allowed by 
BBN |T3. However, the constraints are quite stringent 
for the case of the bino NLSP, and result in an upper 
bound to 3/2 < 10 MeV for M NLS p = 150 GeV. At a fixed 
value of the NLSP relic abundance, the BBN constraints 
give an upper bound on the decay timescale; hence, the 

5 /2 JT^l 

above limit on m 3 / 2 scales as M N ' LSP , see Eq. (JT5J). The 
BBN limit on m 3 / 2 evolves as follows with respect to 
the bino annihilation cross-section; for reference, we re- 
call that the cross-section used in the calculations re- 
ported in Fig. is (onlspv) — 00 (-Mnlsp/100 GeV) -2 
with cr = 10- 9 GeV- 2 and M NL sp = 150 GeV. If cr is 
decreased by a factor 10 to 100, the upper bound on m 3 / 2 
shifts to ~ 100 MeV; if, conversely, the cross-section is in- 
creased by a factor 10 to 100, the upper bound on m 3 / 2 
shifts to ~ 3 MeV. As mentioned in the previous section, 
we have implicitly assumed a branching ratio to hadronic 
decay -Bhad = 0.15; if the bino is nearly degenerate in 
mass with Z, the hadronic decay mode is suppressed, 
with a value possibly as small as -Bhad ~ 10~ 3 . In this 
case, for <7o chosen as above, the bound on m.3/2 would 
be ~ 100 MeV, increasing to ~ 10 GeV if 00 is increased 
by a factor 100, and remaining constant if <7o is decreased 
by a factor as large as 100. Overall the BBN constraints 
result in a bound 7713/2 ^ 10 MeV — > 1 GeV depending 
on the bino mass, annihilation cross-section and hadronic 
branching ratio. 

Figure illustrates the so-called light gravitino prob- 
lem: if TO3/2 ^ 1 keV, gravitinos overdose the Universe 
and/or disrupt BBN unless the reheating temperature 
is low, T RH - min [l0 8 TO 3/2 (M 3 /lTcVr 2 , M x /10], and 
gravitinos are not too heavy, m 3 / 2 ^ 10 MeV — 1 GeV. 
Although it is not impossible to achieve such small re- 
heating temperatures, either by low-scale inflation or a 
late phase of thermal inflation, it is not particularly at- 
tractive as it puts further non-trivial requirements on the 
model and may pose problems for a successful genesis of 
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FIG. 1: Contours of Q3/2 in the plane Trh — rn 3 / 2 f° r "fast" messenger decay, as discussed in the text. White corresponds to 
^3/2 > 1; increasingly lighter shaded areas (blue, red and green respectively) from left to right (in the left panel) correspond 
to hot, warm or cold dark matter respectively with 0.01 < O3/2 < 1. The lightest zone (yellow) indicates 03/2 < 0.01. The 
area shaded (in red) by lines oriented NE-SW is the zone excluded by BBN constraints on NLSP decay. The zone shaded with 
horizontal lines is theoretically excluded as it corresponds to Fs > F. The NLSP is bino-like in the lower panels, and stau-likc 
in the upper panels. The messenger mass scale is Mx = 10 5 GeV in the left panels and Mx = 10 10 GeV in the right panels. 



baryon or lepton asymmetry. Moreover it is necessary to 
tune the reheating temperature to the gravitino and/or 
messenger mass, eventhough these quantities derive from 
sectors of the theory that are a priori unrelated. 

The region TO3/2 ^ 1 keV is devoid of constraints on 
Trh since the gravitino is so light that even at ther- 
mal equilibrium, it cannot overdose the Universe. How- 
ever such light gravitinos make up dark matter that is 
too warm to reproduce existing data on the large scale 
structures. Heavier gravitinos rn.3/2 ^ lOMeV are ex- 
cluded by big-bang nucleosynthesis constraints if the 
NLSP is a bino. Since NLSPs decay at time tnlsp ~ 
6 x 10 4 sec (A/ NL sp/100GeV)- 5 (m 3/ 2/lGeV) 2 , the heav- 
ier the gravitino the later and the more constrained the 
decay. Figure ^ differs from those shown in Refs. 
because we have included constraints from GMSB phe- 
nomenology (namely k < 1) as well as updated con- 
straints from BBN and structure formation, and a more 
accurate calculation of the gravitino relic abundance. 
Overall, one finds that the range of allowed m 3 / 2 is 
severely restricted when compared to previous studies. 



In particular, the present conclusion is at variance with 
Ref. [3 which argued that for a sufficiently small messen- 
ger mass scale Mx ~ 10 5 GeV and sufficiently large grav- 
itino mass m 3 /2 ^ 2 GeV, it is possible to find solutions to 
the gravitino problem for arbitrarily high reheating tem- 
peratures. The discrepancy with [6j is tied to the neglect 
in that study of the SUGRA induced MSSM particles 
contribution to the gravitino abundance at large reheat- 
ing temperatures, as discussed in Section III. A., as well 
as of the big-bang nucleosynthesis constraints on NLSP 
decay. As can be seen in the left panels of Fig.Q] the up- 
per bound on Trh does indeed shift upwards, albeit for 
larger m 3 / 2 than expected in [6j, i.e. tti 3 / 2 > 10 GeV for 
Mx ~ 10 5 GeV. This m 3 / 2 region is however forbidden 
by BBN constraints on energy injection, for both stau 
and bino NLSP. If the gluino mass scale is smaller than 
the fiducial value of 1 TeV, say M 3 = 300 GeV, this region 
of parameter space where the bound on Trh is relaxed, 
moves to smaller m 3 /2, i.e. m 3 /2 ^ 0.4—1 GeV. However, 
the BBN constraints also move to smaller m 3 / 2 since the 
NLSP mass is reduced as the gluino mass. Finally, since 
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the MSSM particle contribution to f2 3 / 2 contains a de- 
pendence on m 3 / 2 and Tpjj, we find that Trh is always 
bounded from above due to the BBN bound on m 3 / 2 ; 
for instance, for M3 = 300 GeV and for a stau NLSP, 
the maximal reheating temperature where CI3/2 < 1, is 
~ 10 9 GeV. In the end it turns out that some fine-tuning 
betwen M3, Mnlsp and m 3 / 2 is required to find a region 
in which Trh can become as large as ~ 10 9 GeV. 

There exist other potential renormalizable interaction 
terms that violate messenger number by one unit. (Such 
operators lead to 1/mpj suppressed proton decay [Isj . 
due to gauge coupling of the messenger fields.) In the fol- 
lowing we assume for definitencss that messengers come 
in complete representations of SU(5), in particular, as 
5m + 5m (or 10a/ + 10a/), while the visible sector su- 
perfields are denoted by 5 f + Wf, and the Higgses sit in 
one pair of 5// + 5/r and one 24// supermultiplets. Gauge 
symmetry limits those interaction terms between messen- 
gers and visible sectors superfields in the superpotential 
to the following: 



W lcn D { 5m5f,h10f , 5 A f 10^10^ , 5 A /5i? ff 24ff , 
5a/5// 24// , 10a/5//5// , Wm5h,f5h,f , 
10a/10 f 5 h , 10 f T0 m 24 h }. (23) 

The interaction terms considered in Ref. 01 are con- 
tained in 5a/5#10 f and 10a/10/?5//, and lead to fast 
messenger decay. Slow decay may, in principle, result 
from operators in Eq. (|23|l which involve particles with 
GUT scale masses. However, renormalizable couplings 
to a 24// must be excluded as they lead to Mgut mass 
mixings between messengers and visible sectors particles, 
hence they would spoil the phenomenology at the elec- 
troweak scale. The colored triplet GUT Higgs in 5h and 
5h carry masses of order Mgut (but no vev), but they 
do not couple to the lightest messenger, which is either 
the 5 th component of a 5m or the (4, 5) component of a 
10a/. Hence the colored Higgses do not lead to suppres- 
sion of the decay width. In Ref. Q it was proposed that 
delayed messenger decay could occur if the decay of the 
lightest messenger was suppressed by the mediation of a 
particle of mass ~ 10 12 GeV in a renormalizable interac- 
tion. Our present discussion shows that this model is not 
natural in the sense that it requires a new particle with 
both the required mass ~ 10 12 GeV and gauge charges 
such that the required renormalizable interaction could 
occur. In the above list of possible interactions terms, 
such coupling does not appear for the minimal content 
of 517(5). 

One may also argue that the required fine-tuning to 
avoid flavor changing neutral currents (in particular for 
light messengers) may actually indicate the absence of 
those renormalizable interactions, and that messenger 
number violation occurs via further suppressed interac- 
tions. Such couplings will be discussed further below. 



2. Renormalizable couplings in the Kdhler function 

Fujii & Yanagida [j| have proposed messenger-matter 
mixing due to a correction in the superpotential SW ~ 
((IU)/m 2 31 )5A/5F, where (W) ~ m 3 / 2 r7i|j. In the frame- 
work of supergravity, a possible origin for such a superpo- 
tential term can be highlighted by adding to the minimal 
Kahler potential Kq — J^. (i running over all su- 

perfields <i>i) a non- minimal part SK given by, 



SK = 5 a/5/? + h.c. 



(24) 



SK is allowed by gauge symmetries (and possibly by an 
R-symmetry as well, for conveniently chosen R charges). 
Then, making use of the usual invariance of the super- 
gravity Lagrangian under Kahler transformations K — > 
K + F($) + followed by superpotential W -> 

e~ F W (and super- Weyl) scalings, the above SW is ob- 
tained for F($>) = —SK to the lowest order in l/mpj, 
provided that a constant is added to the superpotential to 
fine-tune the cosmological constant after supersymmctry 



breaking (whence (W) 
the lightest messenger 



: m 3 /2 m p\) ■ I n our notations, 
is a linear combination of the 



lightest scalar components of 5a/ and 5 M , see Section 2.1. 
The mixing between 5m and 5p generated from SW ~ 
™3/25a/5f thus leads to the decay of X into a lepton and 
a gaugino with width Tx^ix — (g 2 /16n) my 2 /Mx 

We should stress here that, starting as we do from SK 
rather than SW of one expects further contributions 
to the messenger decay or annihilation, with possibly im- 
portant effects on the final gravitino abundance. Indeed, 
other contributions to the decay into visible sector parti- 
cles originate from the supergravity scalar potential pflj 



V B 



K ij W 



Ki 



rn 



Pi 



Wi)[W 



K 



j_ 

l-py 



w;. 



3WW* 



pi 



3* 



(25) 



K %3 



where, i,j* label the full set of scalar fields 
denotes the inverse of the matrix dK/dfadfij, , and Wi — 
dW/dfa, Wp = dW*/d(j)*,. From K D K + SK and 
taking for illustration the case of 5 + 5 messengers with 
W D 55a/5a/ + ?/5_f5//10f + (W) one finds the leading 
contributions to the potential which induce the decay of 
the lightest messenger, 



V B D m 3 f 2 S*5 F 5 M + ym 3 /25 M 5* H lO* F - 2ml /2 5 M 5 F 



m 



Pi 



5a/5f 



\S\ Z (5 M 5* M + 5a/5a/) 



Altogether this is very reminiscent of the Giudice-Masiero mech- 
anism which provides a solution to the so-called /^-problem |44| . 
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+5m$m\ 



y 2 \w y 



+h.c. 



2/5 M 5 



h 10* f )(S5 m 5m 



y5 F 5 H 10 F ] 



keep those contributions under control, eventhough they 
originate from the gravitational non-renormalizable sec- 
tor. The companion operators involving the messenger 
fcrmions are contained in the 0(l/m pl ) part of 



(26) 



where we have neglected terms suppressed by higher pow- 
ers of l/ra pl or of order m 3 /2/ m pi an( A smaller. Other 
operators which do not induce messenger decay into stan- 
dard model particles are not shown, being irrelevant for 
the present discussion. After the spurion scalar field S 
has developed a supersymmetric vev, the first term of 
Eq. H26|) leads to the bilinear operator m^^Mx^F^M 
contributing to the decay considered in [8| (actually a 
similar mixing between the fermionic partners is also gen- 
erated, see below). Note that the second operator in 
Eq. (|2(j|) can mediate an equally efficient decay of the 
lightest messenger to a (colorless) Higgs and a slepton 
r ' x ~>Hf — (y 2 /32ir) m 2 / 2 /Mx, provided that the Yukawa 
coupling is large (i.e. if SK involves the third gener- 
ation). A similar decay to a Higgsino and a standard 
model fermion occurs as well, triggered by the MSSM 
Yukawa vertex and m- i j- 1 Mx5F5 M . In any case, these 
new contributions do not lead to a significant change in 
the analysis of the gravitino relic density, being of the 
same order as Tx^ix- 

Other contributions in Eq. I|26|) can potentially lead 
to important modifications when loop effects are con- 
sidered. This is due on one hand to the supersymmetry 
preserving non-renormalizable operators of the form 
5 m ■ ■ -4>4>* ■ ■ -5 f / fnp\, with <ft denoting an arbritrary field, 
and on the other hand to renormalizable operators which 
induce hard supersymmetry-breaking after cancellation 
of the cosmological constant. (In connection with the 
latter operators, the presence of the spurion, a visible 
sector singlet, could destabilize the hierarchy of the 
messenger and/or electroweak scales |45|.) A thorough 
discussion of these issues which are relevant for the 
theoretical consistency of the effective supergravity 
model is out of the scope of the present paper. Hereafter 
we give only a partial and brief discussion of the two 
types of operators. 

After the vev shift S — » S + Mx , the term of order 
in Eq. (|26[) gives the operators 

Ml 



'PI 



These operators generate potentially very large cor- 
rections to the 5m — 5f mixing, through one- loop 
contributions of the 5m and 5m, which are of the order 



of A^ 



where A is the cut-off scale of the physics 



underlying the effective supergravity Lagrangian. Even 
more, the operators containing SS* instead of M x 
lead to more dangerous corrections of order A 4 due to 
two- loop diagrams. The same diagrams can also induce 
direct decays of the lightest messenger into MSSM 
particles, such as X — > j(Z) + I. It is thus important 
to assess the supersymmetric cancellations which would 



D e 



K/2r, 



-V % D Wr R X 3 L + h.c. 



(27) 



where ViDjW = Wy + (K^/m^W + (Ki/m^DjW + 
(Kj/rr&JDiW - (K^/m^W - T%D k W 



^ R ipL S {5 M 5 F + h.c.) 



(see |29|). They read 



£f 3 



2rn P1 



+1p R 5 M TpL^F S 



+ h.c. (28) 



where -0 denotes the Dirac field which combines the two 
fermionic components of the 5 + 5 messenger superfields. 
[The order of occurrence of the fields in Eq. i|28|) in- 
dicates how they are combined into SU(5) invariants 
.] In the supersymmetric limit, the one-loop contri- 
butions to the 5m — 5f mixing induced by — Lp with 
S — > (S) cancel exactly the ones from the scalar loops 
discussed above. After SUSY breaking through the vev 
of the auxiliary field Fs, no dependence on the ultra- 
violet cut-off A is reintroduced. In particular, even log A 
terms cancel out (though they would not have altered 
the size of the mixing) yielding a correction of order 
-Fj/mpj in the limit Fs -C M x , and of order M^-/m P1 
in the limit Fs — M x , which remains negligible when 
compared to the tree-level mixing magnitude m 3 / 2 Mx, 
Eq. (|26(l , in the parameter space region relevant for grav- 
itino dark matter. 6 Finally, as mentioned before, some 
hard SUSY breaking operators are generated in Eq. I|27|) 

from m 3 /2e- ft: / 2mpl (^XRXL^^i+^- c -); leading to bilinear 
matter fermion mixing between the messengers and the 
MSSM particles, m3/a(/ R ipL + h.c). These can poten- 
tially lead to quadratic divergences which would desta- 
bilise the scale of the mixing 5m — 5f between the light- 
est messenger and the MSSM scalar particles. How- 
ever, leading one-loop (tadpole) effects with one mass 
insertion occur as corrections to 5* M ~ 5p and turn out 
to be at worse 0((A 2 /mp 1 )m 3 /2Mx), thus harmless for 
a cut-off of order ~ mpj. (Note also that the matter 
fermion kinetic terms induce fermionic tadpole contribu- 
tions which cancel among themselves due to the deriva- 
tive couplings.) To summarize, in the sector of the su- 
pergravity Lagrangian not involving the gravitino the SK 
piece of the Kahler potential leads essentially to the same 



6 We checked also for cancellations to two-loop order considering 
subclasses of Feynman diagrams which involve the spurion scalar 
and fermion virtual contributions. A detailed study is outside the 
scope of the present paper. 
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features as discussed in [8j, at least to leading order in 
l/mp| and up to one-loop. One exception is the messen- 
ger decay induced by first term in Eq. I|26J1 which we will 
consider later on. 

The relic abundance for the gravitino can be calcu- 
lated using the techniques developed in Section III and 
the results are shown in Fig. [2] This figure uses the color 
shading as in Fig.^but is plotted in the plane m 3 / 2 — A/x 
instead of Irh, which was taken to be Trh = 10 12 GeV. 
The plot shown in the left panel assumes that the spu- 
rion S is much heavier than the lightest messenger X, 
in which case annihilation XX* — > GG takes place 
with the cross-section given in Eq. 10. As discussed 
in Section II. B, this cross-section violates unitarity for 
Mx > 10 7 GeV(TO 3 / 2 /l keV) 2 / 3 , i.e. in the region above 
the dashed thick line in the left panel of Fig. [21 There is 
no solution in this case for gravitino dark matter, at vari- 
ance with the conclusions of Ref. § . The main reason is 
that annihilation into goldstinos has not been accounted 
for in Ref. Q, yet the solution for gravitino dark mat- 
ter proposed by these authors lies in the region in which 
unitarity is violated. The annihilation cross-section at its 
unitarity bound is much larger than that used in Ref. @ 
for annihilation through gauge interactions, hence the 
messenger relic abundance and the amount of gravitino 
dilution are correspondingly smaller. At the very least, 
since one cannot predict the cross-section in this region 
where multi-goldstino production violates unitarity, one 
can conclude that the results for the scenario of Ref. Q 
in this region are model dependent in that they require 
contributions from the hidden sector to bring down multi- 
goldstino production to a negligible level. 

In the right panel of Fig. [21 it is assumed that S is 
much lighter than X. Annihilation into goldstinos scales 
with the effective Yukawa coupling F$/F of the messen- 
ger to the goldstino component of the fermionic part- 
ner of S, as in Eq. ©. This annihilation channel thus 
contributes only in the region of direct GMSB scenarios 
where k = Fs/F ~ 1, which happens to be that where 
gravitino dark matter can be found for the mixing term 
5jv/5f proposed in Ref. The lightest messenger can 
also annihilate into a pair of spurions, as discussed in 
Section II. B. Furthermore, one must also consider the 
possible decay of the lightest messenger into one visi- 
ble sector sneutrino and two goldstinos, which is induced 
by the above mixing term and XX* GG four-point ver- 
tices similar to the ones discussed after Eq. ©. For in- 
stance, the latter vertex is induced by terms of the form 
-e K / 2m ^W^ L a^W R /m 2 !, 29], using Eqs. J2H© and 
W -> (W). 

One expects this decay width to scale as r ' X ->ogg ~ 
0.1(Fs /F) 4 Tx^i\, with the prefactor of 0.1 account- 
ing for the enlargement of phase space. This decay 
channel produces highly relativistic gravitinos which mix 
with the "cold" gravitinos produced by other channels, 
see Section III. A, resulting in mixed dark matter. The 
small area around ro 3 / 2 ~ 10 — 100 keV is shown in red 
(medium shading) in Fig. [21 indicating that the hot grav- 



itinos contribute to more than 10% of the gravitino en- 
ergy density and that the averaged velocity corresponds 
to warm dark matter. Finally one must also include 
decay into a spurion and a sneutrino induced by the 
three scalar coupling, first term in Eq. I|26l) . with width 
T x ^so =s (l/167r)m 2 /2 /Af x . 

The inclusion of these various annihilation and decay 
channels modify the region of parameter space where one 
can find gravitino dark matter with respect to the conclu- 
sions of Ref. 0. The various effects add up in pushing 
this region to higher values of Mx- In effect, the in- 
creased lightest messenger annihilation cross-section de- 
creases their relic abundance hence the amount of en- 
tropy production; similarly, the increased lightest mes- 
senger decay width increases their decay temperature 
hence also decreases the entropy production. These ef- 
fects can be compensated, at a given value of 1713/2, by 
increasing Mx since (crxx'^ v) cx 1/M X and Tx^ oc 
l/M x . 

Note that Fig. [21 assumes that the NLSP is a bino; as 
seen previously, the BBN bounds would be relaxed is the 
NLSP turned out to be a stau. 

Finally, the influence of the choice of the post- 
inflationary reheating temperature (taken as Trh = 
10 12 GeV here) is as follows. For this chosen value of 
Trh, gravitinos (i.e. their 1/2 component) abundances 
are brought to equilibrium for all to 3 / 2 ^ 100 MeV. For 
smaller values of Trh, the plot would thus look similar 
in the region where Mx 7rh, unless Trh is that low 
Trh < 10 12 GeV(m 3/2 /0.1GeV) 2 that gravitinos are ini- 
tially not in chemical equilibrium. For Mx ^ Trh, mes- 
sengers are not produced at the post-inflationary reheat- 
ing, there is no entropy production (gravitino dilution) 
and the gravitino abundance in the plane m 3 / 2 — 7rh 
then resembles that shown in the right panels of Fig. ^ 

If the messengers sit in 10 + 10 representations of 
SU(5), a similar mixing term IOmIOf can be induced, 
and leads to similar effects. Note that the mixing 5m5f 
and 10a/ 10 f are constrained by R— symmetry charge 
assignments. If a non-holomorphic mixing 5 m 5f (or 
10* AI 10p) is allowed, the situation is quite different. In 
effect, the lightest messenger now mixes in the kinetic 
terms with the sneutrino and its decay is no longer sup- 
pressed by to 3 / 2 /Mx- As a consequence, decay is much 
faster, T ~ 0{M X ), and no entropy production occurs; 
the situation is then as shown in Fig. ^ 



B. Non-renormalizable couplings 

1. Superpotential interactions 

It is possible that interactions of the type Eq. (|2*2^l are 
forbidden and that messengers may only decay via non- 
renormalizable operators in W, as discussed previously. 
Operators which may arise due to unknown Planck-scale 
physics, which respect the SU (5) gauge symmetry and 
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FIG. 2: Contours of in the plane Mx — rn 3 / 2 f° r one pair of messengers sitting in 5 + 5 representations; the lightest 

messenger X carries the gauge charges of a vl, and decays into a slepton and a gaugino via the mixing term 5a/ 5f in the 
Kahler function. Color shading is as in Fig. the NLSP is assumed to be a bino in both_panels. In the left panel, the spurion 
S is assumed to be much heavier than the lightest messenger and annihilation XX* — » GG proceeds with the cross-section of 
Eq. |7] the dashed line indicates the limit above which this cross-section violates the unitarity bound, see Section [H] In the 
right panel, the spurion S is assumed much lighter than X, and annihilation into goldstinos is given by Eq. |S| See text for 
further details. 



that violate the messenger number by one unit are given 
to leading order in 1/mpi, by: 

Wnon-ren D { 5m10 F 10 F 10 F , 5 M 5h5h,f5h,F , 

5 m 5h5h5h,f , 5m5h5hWf , 
5jvf5ijl0i?24ij , 5m5h,f24:h24:H , 
5M5ff24ff24fj , 10f10m5h5h,f , 
10m5h,f&h,f5h,f j 10m5h5h5h , 
10 m 5^f10 f 10 f , 10 M 10 F 5 ff 24 H , 
10 M 5_ff 5 H 2A H , 10m5h,f5h,f24h > 
10m10f24# 24h , 5m5a/5a/5f , 
10 m T0m10m10jt} (29) 

All terms in Eq. 129|) which involve couplings of one 
lightest messenger to 5f or 10 p but not Higgses lead 
to decay into three-body final states with decay width 
~ 10 _4 Mj^/m P1 . It is easy to see, using Eq. (|TU|) that en- 
tropy production is not sufficient to dilute the gravitinos 
to the required abundance for a high post-inflationary re- 
heating temperature T RH > 10 8 GeV (and T RH > M x ). 
Admittedly this is a drawback of the present scenario 
since those terms are the most generic. 

Let us now consider the terms involving couplings 
to Higgses. For terms involving one 24# acquiring 
vevs of order ~ Mqut, say A<5>i<f > 2 i?24/"^Pi, the La- 
grangian contains the effective Yukawa interaction ~ 
(Mqut /mpi)XxiX2 horn the fermionic part of the La- 
grangian contained in Eq. I|27|l , with \i an d Xi the Weyl 



spinors of <f>! and <f> 2 - This Yukawa interaction between 
the lightest messenger and two fermions with effective 
coupling constant ~ 1CP 3 leads to fast decay if \i an d 
X2 have electroweak scale masses. Other terms in the 
Lagrangian lead to similar or only somewhat smaller par- 
tial widths. Consequently the conclusions of the previous 
section with regards to gravitino dark matter with " fast" 
decaying messengers apply. From Eq. (|29|l one can check 
that all possible above combinations involving one 24# 
contain a coupling of X to particles with electroweak 
masses except 10Ml0i?5#24# . However, even for the 
latter term, the scalar potential contains the interaction 
10 M 10F24 i? ylO F 10 F /rn P i generated by \dW/d5 H \ 2 , 
with y the third family Yukawa coupling. This interac- 
tion gives a decay width T ~ 10~ 5 ?; 2 (A/gut/topi) 2 ^x, 
which is too large to allow solutions for gravitino dark 
matter. 

Terms involving two 24 should be excluded as they 
lead to unacceptable mass mixings between messengers 
and visible sectors superfields. 

Consider now terms of the form Y$i$2 < I ) 3/'Tipi con- 
taining at least one 5# but no 24#. It can be checked 
that all terms of this form in Eq. I|29|l contain couplings 
of X to particles with electroweak masses, hence lead 
to fast decay as above, except for 5jv/5.jf5#10 F /mpi, 
10m5jj5k5h/"Ipi and 1Q m ^h,f^h,f^>h,f /mp\. The 
first of these terms, when written for the lightest mes- 
senger, contains at least one colored Higgs, say Hi. How- 
ever, the scalar potential term \dW/dbHi \ 2 generates 
here as well an interaction between the lightest mes- 
senger with 4 particles of electroweak masses, leading 
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to decay width T ~ 10~ 5 y 2 M|-/mpj, still too large for 
dark matter solutions. The term 10M5#5#5#/mpi cou- 
ples X to the three colored Higgses, hence its decay 
is too highly suppressed both by the GUT scale and 
phase space, T ~ 10 _12 (Mx/MGUT) 8 Af|-/7n P1 , an d can- 
not lead to decay before BBN for M x < 10 14 GeV. Note 
that the big-bang nucleosynthesis constraints on NLSP 
decay can be translated into an extreme upper bound 
M x ^10 12 GeV (see Fig. and H3). Finally for the 
term 10^/5^.^5^ similar conclusions apply 

if three Higgses are involved. For couplings involving two 
colored Higgses, decay is also too highly suppressed. In 
effect, the lightest messenger can then decay into 5-body 
final state with mediation by a GUT mass Higgs, leading 
to T ~ lCrV (M x /Mgut) 4 M x /m P1 . It cannot decay 
before BBN if M x < 10 12 GeV, as can be checked using 
Eq. I|16fl . Finally, if only one Higgs is involved in the 
coupling, the messenger can decay into four particle fi- 
nal states with decay width T ~ \Q~~ 5 y 2 M x /m 2 ^ and yet 
provide no solution for gravitino dark matter. 

The last two operators of Eq. I|29|l together with the 
first term of Eq induce at one-loop order a mixing 
between the lightest messenger and the MSSM matter 
fields. This mixing is non-vanishing only after super- 
symmetry breaking and is found to be of magnitude 
~ (V2Mx/m P i) F s log(A/Afx). As discussed in Sec- 
tion IV. A. 2, A is identified with a physical cut-off of order 
topi, however, in contrast with the results of that section, 
there is here no (accidental) cancellation of the cut-off de- 
pendence and no suppression by the gravitino mass. The 
mixing is so large that decay is not accompanied by en- 
tropy production, and consequently there is no solution 
for gravitino dark matter. 

In summary, non-renormalizable interaction terms in 
the superpotential for the minimal content of GMSB sce- 
narios in SU(5) grand unification do not allow for nat- 
ural solutions leading to gravitino dark matter (unless 
the post-inflationary reheating temperature is tuned as 
before) . 



2. Non-renormalizable interactions in the K'dhler junction 

Interaction terms in the Kahler potential K can either 
be holomorphic or not in the superfields, leading to dif- 
ferent phenomenologies, which we explore in turn. To 
leading order in 1/mpi, non-renormalizable holomorphic 
operators have the same form as those shown in Eq. , 
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(30) 



We write generically these operators as X4>\<f>2/m?\. Let 
us assume for the moment that both <f>\ and 4>2 have 
masses of order of the electroweak scale. 

The Kahler U(l) connection Kjd^4>j — h.c, with Kj = 
dK/d(j)j induces the following coupling to the fermionic 
components -0 1 of all the chiral superfields of the model 



c f D tt^^lI^l Im d^Xfafa), (31) 
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which leads to a partial decay width cx M x jm\ v From 
Eq. I|27(l an effective Yukawa coupling is generated, 

d 2 K (WW a m 3/2 _i 



nn — —^nth 



X<P R r L , (32) 



which leads to a highly suppressed partial decay width 
~ ( m 3/2/ m Pi) 2 M x /mp l . Finally, couplings to goldstinos 
are generated notably by the gravitino mass term and 
gravitino kinetic terms: 
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_ (33) 
with the replacement ^ M — > 9 A1 G/m 3 / 2 . 

All terms lead to extremely slow decay if <f>i and <p2 
have electroweak scale masses, and must be forbidden in 
order for the lightest messenger not to decay after BBN. 
However if one 24# is present, the replacement of this 
field by its vev in the Kahler function shows that one re- 
covers a mixing term as proposed in Ref. [8| albeit with 
effective coupling Mqut/topi- This mixing term then 
leads to the same decay widths into one sfermion and 
one gaugino, or one sfermion and two goldstinos, as dis- 
cussed before, albeit decreased by (Mgut/^pi) 2 ~ 10 -5 . 
The consequences for gravitino dark matter are shown in 
Fig- El As discussed before, we should consider the cases 
where the spurion is heavier or lighter than the light- 
est messenger. The left panel shows the case where the 
spurion is heavier than the lightest messenger and the 
annihilation cross-section into goldstinos increases with 
increasing Mx to saturate at the unitarity bound above 
the dashed line. Only a small portion of parameter space 
allows for gravitino dark matter in this case; it is ac- 
tually an amputated part of the solution shown in the 
right panel, see below. Above the dashed line, the re- 
sults are highly uncertain since multi-goldstino is not 
well controled. Moreover, the lightest messenger can de- 
cay into a sneutrino and a pair of goldstinos due to the 
above mixing, with a decay width which is expected to 
scale as r x ^gg k { F l M x I fA ) t x^ix- The prefactor 
(FjM^/F 4 ) denotes the effective coupling of XX* to a 
pair of goldstinos in the heavy spurion limit. In the re- 
gion in which unitarity is violated, it has been assumed 
that this effective coupling saturates at the value reached 
[~ 0(1)] when the annihilation cross-section reaches the 
unitarity bound. It is then comparable to the decay 
width into a lepton and a gaugino and produces highly 
relativistic gravitinos. The energy density contained in 
these gravitinos exceeds the BBN bounds on additional 
relativistic degrees of freedom so that most of this region 
is excluded, as indicated by the NW-SE oriented dashed 
lines in Fig. O The SW-NE oriented dashed lines exclude 
the part of parameter space at small m 3 / 2 in which the 
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FIG. 3: Contours of ^3/2 in the plane Mx-m 3| / 2 for one pair of messengers sitting in 5+5 representations; the lightest messenger 
X decays into two goldstinos and one sfermion or one sfermion and one gaugino via the mixing term ~ X0i(Mgut/»tipi) in 
the Kahler function. Same color code and conventions as in Figs. [5] The red NE-SW dashed area in the left part of the 
figure results from big-bang nucleosynthesis constraints on lightest messenger decay and the NW-SE area is forbidden by the 
contribution of gravitinos to the energy density at BBN; see text for further details. 



lightest messenger decay occurs so late that it is forbid- 
den by BBN constraints on energy injection. If the NLSP 
were a stau, this region would still be forbidden but the 
constraints at large 772,3/2 would be relaxed. 

In the right panel of Fig. |3| the spurion is assumed to 
be lighter than the lightest messenger and annihilation 
into goldstinos is less effective. One finds a solution for 
gravitino dark matter in a large part of parameter space, 
m 3 / 2 ~ lOkeV — ► lMeV, for scenarios of indirect gauge 
mediation, i.e. k <C 1. This solution is the same as that 
discussed in Section IV. A. 2, albeit shifted to smaller val- 
ues of Mx\ this can be understood from the fact that 
for a same relic abundance of X, entropy production is 
larger in the present case since the decay width is further 
suppressed. Hence one can tolerate a smaller relic abun- 
dance, or equivalently a higher annihilation cross-section, 
i.e. a smaller Mx- 

Finally consider now non-holomorphic non- 
renormalizable couplings between X and visible 
sector particles in K. Such couplings can take the form: 

K D 7^r{ 5m5 h ,f10 f , 5 M 5]yl0jr , 5^10 F 10jr , 

5 m 5h,f24h , 10^5^5^, 10 M ^ H F 5 H , 

1o\.j5h,f5h.f , 1o\j10f5h , 10m10f5jj : f , 
10t, f 10 F 24 ff , 10 M 10 f F 24 ff , + h.c.} (34) 

As before, we write this coupling as X* fafa / 'mpi and 
assume for the moment that fa and fa carry electroweak 
scale masses. Then decay into fa an d fa with width 



r ~ 10~ 2 M| r /mp 1 occurs via the mixing of kinetic terms 
between fa and X or between fa and X. As seen be- 
fore, such a decay width does not lead to solutions for 
gravitino dark matter as entropy production is not sig- 
nificant. Inspection of Eq. (|34fl reveals that all terms fall 
in the above category except those involving one 24# as 

well as w\j10f5h and 10j\/10f5^. 

The latter two terms necessarily contain one colored 
Higgs with GUT mass but no vev, which we assume to 
be fa. Then the scalar potential term involving the 
inverse of the Kahler metric g contains a coupling of 
X to visible sector particles: g* 2 M * WD M ,W* D 
(fa/mpi)yMx lOjrlOirX. This leads to decay into three- 
body final state with width T ~ 10~ 4 y 2 M|-/mp 1 , again 
too large to yield solutions for gravitino dark matter. 

Finally, if coupling to one 24# occurs, say fa, mass 
mixing of order MqutMx occurs between fa and X and 
leads to one negative mass squared eigenstate; this cou- 
pling must therefore be excluded. 



V. DISCUSSION 

In this section we explore qualitatively other possible 
avenues which may help reconcile gravitino dark matter 
with more generic GMSB scenarios. Indeed, the previ- 
ous discussion has shown that interesting solutions for 
gravitino dark matter and/or the gravitino problem in 
GMSB with SU(5) grand unification with a high post- 
inflationary reheating temperature, can only be found for 
some very specific couplings between the messenger and 
visible sector and in some restricted regions of parame- 
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FIG. 4: Contours of in the plane Mx — W3/2 for one pair of messengers sitting in 5 + 5 representations as in Fig. [5] but 
for decay width Y ~ lO~ w M x /m% l . See text for details. 



ter space. It is furthermore necessary to assume that the 
spurion is much lighter than the lightest messenger so 
that multi-goldstino production remains at a safe level. 
A gravitational decay width V ~ 10~ 3 Af|-/?7jpj, which is 
generic in the sense that it is generated by most allowed 
non-renormalizable messenger-matter interaction terms, 
does not lead to satisfying solutions for gravitino dark 
matter. 

It is instructive to consider the case of decay widths 
with the same scaling but whose prefactor is much 
smaller, T ~ eM x /m P1 with e <C 10~ 3 . Figure 0] shows 
the solution for e = 10~ 10 . Such decay width can be 
achieved by terms of the form W D X&i&2$3H24/m Pl , 
with H24 designing a Higgs with non-zero vev in 24#, 
or by most non-renormalizable couplings to order 1 / mp\ 
discussed in the previous section, provided they are fur- 
ther suppressed by a factor ~ 10~ 7 . In this figure, one 
finds that in the left panel, where S is assumed heavier 
than X, i.e. where multi-goldstino production plays a 
significant role, there is room for gravitino dark matter 
only in a very limited region of parameter space. In this 
area, furthermore, the post-decay reheating temperature 
is quite close to 1 MeV. On the contrary, in the right 
panel one recovers a solution for gravitino dark matter 
with mass TO3/2 ~ 10 keV — > 10 MeV in direct gauge me- 
diated scenarios k < 1. Higher values of the coupling e 
lead to solutions shifted to higher Mx, with the dashed 
region due to late messenger decay shifting downwards. 
Lower values of e lead to solutions shifted to smaller val- 
ues of Mx , but with the excluded dashed region moving 
upwards in Mx- 

Dangerous operators involving GUT-scale vev's may 
also be present. However, global continuous R- 
symmetries are expected to play an important role in sce- 
narios of supersymmetry breaking [46j, and in a generic 



setting they could control the absence of unwanted op- 
erators for properly chosen R-charges 01 ■ Discrete 
Z-symmetries, motivated by the need to improve the 
fine-tuning issues can also play a selective role. 

In particular, the spurion gauge singlet superfield S 
can be present or not in non-renormalizable operators 
in the Kahler or superpotential depending on its R— 
and Z— charges attributions. Such operators, provided 
that they do not destabilize the mass hierarchies 
pUf. would lead to tree-level suppressions of the form 
{(S)/m P i) 2n ~ (Mx/mpi) 2n . However the decay width 
is now too suppressed to yield reheating before BBN ex- 
cept in a very narrow region centered on m 3 / 2 ~ 1 MeV 
and M x ~ 10 10 GeV. 

In contrast, a larger number of seemingly generic so- 
lutions to the gravitino problem/gravitino dark matter 
may be found if one considers S'O(IO) grand unification 
and Ms > Mx- This case is studied analytically in a 
companion paper 26] for Mx ~ 10 6 GeV. In Section^ it 
has been shown that the amount of entropy produced de- 
pends directly on the relic abundance of the lightest mes- 
senger, which in turn depends directly on its nature. In 
50(10), the lightest messenger is a SU(3) x SU(2) x U(l) 
singlet (ur— like), hence its annihilation cross-section is 
suppressed: it may either annihilate through one-loop di- 
agrams or at tree level into goldstinos, and (at tree level) 
through suppressed diagrams of GUT mass bosons ex- 
change. One may estimate the relic abundance of the 
lightest messenger in this case by using the dimensional 
estimate (<Jxx^...v) ~ (a/^Tr) A /M x for one-loop dia- 
grams and the annihilation cross-section into goldstinos 
computed in Section II. B (see also Ref. |2(J). One finds 
that the relic abundance is larger than for the SU(5) case, 
hence the amount of entropy production is expected to 
be correspondingly larger. One then finds that a generic 
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gravitational decay width T ~ 10~ 3 M^-/mp 1 can lead to 
natural solutions for gravitino dark matter in a signifi- 
cant part of parameter space, as shown in Fig. El 

This solution is attractive for several reasons. As seen 
in the left panel of Figs. El gravitino cold dark matter 
with the right relic abundance occurs in a rather large 
part of parameter space, hence it appears "natural" in 
this sense. Moreover the decay width to sparticles as- 
sumed T ~ 10 _3 M|-/TOpj is generic as it is predicted 
by most non-renormalizable operators that violate mes- 
senger number by one unit. In this region of parameter 
space, the factor k = F$/F sa 10~ 5 — 10~ 2 , hence grav- 
itino dark matter would be obtained for indirect gauge 
mediated scenarios; in Ref. j2fj it is argued that this 
case can be naturally incorporated in the simplest indi- 
rect GMSB scenarios |2[. Furthermore, the solution for 
gravitino dark matter occurs in the predictive region in 
which multi-goldstino production satisfies the unitarity 
bound, unlike the solutions seen hitherto for SU(5). Fi- 
nally, Fig. El assumes, as the previous figures, that the 
NLSP is a bino; for a stau NLSP, the BBN constraints 
at large m 3 / 2 would be relaxed, and this would enlarge 
in turn the space of solutions for ^3/2- 

One can show |2q that the amount of entropy pro- 
duction does not hinder successful leptogenesis at high 
reheating temperatures; this is all the more interesting 
as leptogenesis scenarios typically operate in models of 
50(10) grand unification rather than SU(5). Strictly 
speaking, the above scenario requires the spurion field to 
be heavier than the lightest messenger. This issue is how- 
ever model-dependent as was briefly discussed in section 
II. B. For completeness, we illustrate in the left panel of 
Fig- El the opposite configuration where the lightest mes- 
senger annihilation into a pair of spurion fields is con- 
trolled by Eq.JSJ. This annihilation leads to a too low 
messenger relic density for the entropy dilution mechan- 
sim to work. 

Finally we note that the results obtained in the 
present study remain valid when R— parity is violated. 
In effect, in this case the gravitino lifetime is t 3 / 2 ~ 
10 20 sec (TO3/2/I GeV)~ 3 for trilinear R— parity violating 
terms 15] or r 3 / 2 ~ 10 27 sec (ma/2/l GeV)~ 3 for bilin- 
ear R— parity violating terms |14| . Hence the gravitino 
is sufficiently long-lived that it can be considered as sta- 
ble dark matter with respect to the formation of large- 
scale structure. If the gravitino lifetime > 10 27 sec one 
also finds that distortions of the diffuse backgrounds 
due to gravitino decay are evading observational con- 
straints. For trilinear R— parity violating terms, this re- 
quires TO3/2 ^ lOMeV, while for bilinear terms, TO3/2 ^$ 
1 GeV is sufficient. With regards to the NLSP, its de- 
cay can proceed into visible sector particles on a short 
timescale and BBN constraints can be evaded, albeit they 
are replaced with constraints on diffuse background dis- 
tortions. Hence the plots in parameter space would look 
similar to what has been found above. 

VI. CONCLUSIONS 

We have presented an exploratory though detailed in- 
vestigation of relic LSP gravitino abundances in scenar- 



ios of gauge mediated supersymmetry breaking (GMSB). 
This study focuses on the possibility of gravitino dark 
matter and on solving the light gravitino overproduction 
problem for reheating temperatures after inflation that 
are " arbritrarily" high. GMSB scenarios contain inter- 
mediate mass scale 10 4 GeV < Mx ;$ 10 12 GeV messen- 
ger fields which by virtue of their gauge interactions are 
easily produced in the primordial plasma. Cosmology re- 
quires these particles to subsequently decay as they would 
otherwise overdose the Universe (except for a lightest 
messenger with Mx ~ 10 — 30TcV). Flavor-changing 
neutral currents impose somewhat restrictive limits on 
messenger number violating Yukawa interactions, possi- 
bly arguing for such messenger number violation to be 
rather weak. If so, the delayed decay of messengers 
may subsequently dilute any pre-existing gravitino abun- 
dances in accord with cosmological constraints. 

We have thus investigated a fairly complete set of 
renormalizable and non-renormalizable messenger num- 
ber violating operators within supersymmetric unifica- 
tion in 5/7(5) (as well as some within 50(10)) and 
their impact on relic gravitino abundances. Results are 
shown for a variety of operators and imposing relevant 
constraints on NLSP decay and messenger decay from 
BBN, as well as constraints on the "warmness" of grav- 
itino dark matter from the required sucessful formation 
of large-scale structure. With respect to prior, less de- 
tailed, studies 0, IE U El IbH ITsj- we have uncovered a 
number of significant changes, notably the importance of 
messenger-messenger annihilation into two goldstinos in 
part of the Mx - TI3/2 parameter space, which modifies 
the messenger pre-decay freeze-out abundances in SU (5) 
and 5O(10) grand unification. 

In general, we have found that gravitino dark matter 
in 5t/(5) grand unification in scenarios with high post- 
inflationary reheating temperatures Trh is only possible 
for a few specific messenger-matter couplings. Further- 
more we have shown that these models predict gravitino 
dark matter in regions of parameter space in which mes- 
sengers annihilation to goldstinos violates unitarity un- 
less one makes specific assumptions on the mass spectrum 
of GMSB models, and in particular, that the spurion 5 
be much lighter than the lightest messenger. 

In contrast, in 5O(10) grand unification gravitino dark 
matter may be obtained for a variety of generic operators 
and in the predictive region of parameter space where 
multi-goldstino production is under control, as long as 
renormalizable messenger number violating interactions 
in the superpotential are absent [26|. We thus believe 
that gravitino dark matter in GMSB scenarios is a viable 
alternative to neutralino (and gravitino) dark matter in 
supergravity scenarios, and as such deserves further de- 
tailed study. 
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FIG. 5: Contours of in the plane Mx — TI3/2 f° r one pair of messengers sitting in 16 + 16 representations of 5*0(10); 
the lightest messenger X is a singlet under SU(3) x SU{2) x U{\). Its loop-suppressed annihilation cross-section scales as 
(a/4iv) 4 /M x , and it decays into sparticles through non-renormalizable operators with width T ~ 10 -3 Mf/m|j. Color shading 
is as in Figs. [T] H See text for details. 
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